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We present a fundamental solution to an initial value problem for the KdV- 
Whitham system in an explicit integral form. Monotonically decreasing initial 
data with finite number of breaking points are considered. Generating function 
for the commuting flows of the averaged KdV hieararchy producing the analytical 
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solutions to the KdV-Whitham system is constructed. 
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1 Introduction 

It is known [1,2] that evolution of smooth initial distribution 

u(x, 0) = Uo(x) (1) 

according to the Korteweg-de Vries (KdV) equation 

^ i d t u + 6ud x u + e 2 dl xx u = (2) 

is described locally as e — > by g - phase (g - gap) solution which is expressed in terms 
of theta-function of the hyperelliptic Riemann surface of genus g (see Refs. [3,4]) 

29+1 

r : y 2 = n(A-r i )= J R 2g+1 (r,A), (3) 

i=i 

n < r 2 < ... < r 2g+ i. 

This solution has the form 



d 2 , „ ( kx — ujt 



u(x,t) = —2—^ In 9 
ax 



r 1 ,...,r 2g+ i +c(r) . (4) 
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Here the components of the wave number and the frequency vectors k(r) = (k^\ k^j 

ui(r) = (uj^\ ...,00^^ are expressed in terms of the branch points rj by well-known 
"finite-gap" formulas 

k (m) = ffl M ? Jm) = ffl M 2 £ r . + ^ m = 1? ^ (5) 

where the functions ae 4 - m ^(r) are defined uniquely by the normalization conditions 

fdQ^=5 km , k,m = l,...,g (6) 

for the basis holomorphic differentials 

^ M = 1 4 m) (r) ^ d\ , (7) 
k=i yR 2g+ i(r,\) 

where the sign is given by y/T = 1. 

We consider the Riemann surface (3) with cuts along the intervals (— oo, ri], [r 2 , r 3 ], 
...,[r 2g , r 29+ i]. The canonical basis of cycles a\,...,a g ; @i,...,P g is chosen so that the 
cycles a%, a g He on T over the permitted zones [ri, r 2 ], [r 29 _i, r 2s ]. The genus of 
the Riemann surface depends on the initial data uq(x) (globally) and on the position 
on the (x,t) - plane (locally). In the general case the (x,t) - plane is splitted into a 
finite set of regions in every of which the solution is descibed by the formula (Q) with 
certain g [5,6] (Fig. 1). The slow x, t - evolution of the branch points rj parametrizing 
the "uniform" solution (|J) is governed by the Whitham modulation equations [7,8,1] 
and provides the required transition. These equations have a diagonal (Riemann) form 

d t r i + V i {r)d x r i = , i = 1,2, ...,20+1 (8) 

with real and different characteristic speeds (pure hyperbolicity [9]) which will be speci- 
fied below. The existence of global solutions to @ corresponding to the monotonically 
decreasing analytic initial data ([!]) for the KdV equation was proven recently by Tian in 
Ref. [10] where a set of analytic solutions in the form of infinite series were constructed 
with the aid of the generalized hodograph transform [11] and Krichever algebraic - 
geometrical procedure [12]. 

Unlike that of Ref. [10] we propose the direct way to construct the general (nonan- 
lytic) global solutions of (^) in an explicit integral form. This problem has been solved 
to date only in the simplest single-phase case [13-16]. Construction of the global so- 
lutions for the multiphase case is vital both from purely theoretical point of view to 
complete the theory of the multi-gap averaged integrable systems and for the physical 
applications where the complex multiphase nonlinear dissipationless wave behaviour is 
not uncommon (Alfvenic turbulence in the solar wind [17], the Earth's bowshock [18] 
and oth.) 
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2 Generating function for the averaged KdV hier- 
archy 

We start from the potential representation for the characteristic speeds Vj : 

3.mH 

^W = ^y; d^d/drj (9) 

which immediately follows from the "conservation of waves" 

<9 t k(r) + «9^(r)=0, (10) 

considered as a consequence of (|j) (for the single-phase case see Refs. [14,15,19]). 
We notice that the representation @ is valid for any chosen phase number since the 
characteristic speeds Vj do not depend on m. 

Looking for the solution of (§) in the characteristic form 

x-V i (r)t = W i (r) (11) 

one arrives at the linear overdetermined compatible equation system for Wj's [11,20]. 

diWj <),\) ... 

■ i^j. (12) 



Wt - Wj V - Vj 

Construction (pTTj) , ( ]T2| ) is known as the generalized hodograph transformation and 
provides all the hydrodynamic symmetries of i.e. equations of the form 

d T n + Wi(r)d x n = , (13) 

which commute with @, i.e. dt^i = d T tVi. Here r is "new" independent time. Some 
of the solutions of (0) (the uniform ones) correspond to the averaged hierarchy of the 
KdV equation [8,12,20,21]. The system (|T2| ) can be locally parametrized by 2g + 1 
functions of one variable [20]. We remark that this corresponds to the number of 
arbitrary functions r\, (x, 0) in the initial value problem for the Whitham system (Q). 

Now we shall construct a generating function Wj (r; A) for the uniform commuting 
flows to (H) which will give the characteristic speeds of the averaged KdV hierarchy 
as the coefficients of the expansion of Wj(r; A) in powers of A -1 as A — » oo. The two 
first terms of this series are obvious and correspond to the zeroth (linear equation of 
translation) and the first (system (§)) terms of the averaged KdV hierarchy: 

Wf=l, Wj = Vj . (14) 

So we shall seek Wj(r; A) in the form predicted by (|^), (|T2"|): 

Wj(r;\) = 9 ^ ( ™(r;A) m = l (15) 
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where 

w< m >(r;A) = JfcW+£?^(»0 + ... (16) 

A 

Such a generating function for the single-phase case was constructed by Pavlov in [22] 
and has the form ^ 

^ s .p.O; A) = . (17) 



With (|15|), (|Tq), (|TyD in mind we propose the following project: 

S W (r;A) = A 3 / 2 E4 m) , = = A 3/2 ^- , (18) 
fc =i ^R 2g+ i{r, A) ^ 

which has two first terms of expansion as A — > oo coinciding with fll6|) and turns into 
([ITD in the case g — 1. 

Let us check that commuting flows fll5|) with the frequences (|I8|) identically satisfy 
the linear equation system (|12j). 

With this aim in view we consider the following combinations occuring in the left 
part of (0): 



.^ M / " (A-r J )(A-r,) v / J R 29+1 (r,A) 



and 

A 9 + 6^ ) 1 A 9 " 1 + ... + ^ m) 



/2 /^ ^)\_ (m) 



S^H J (A - r,)(A - rO^iCr.A) 



(20) 



Here pjj^(r), (r) are some certain (analytic) functions. The coefficients o^(r) 



6jt (r) can be uniquely defined by the normalization conditions, namely, integrating 
flig|), (pp|) over the a-cycles we arrive taking into account (El), (18) at two sets of 
algebraic equations 

f A' + 4-iA»- + .- + ^> a 0< 

«* (A-r j )(A-r i ) v /i?2« 7 +i(r,A) 

and (21) 
A' + &£ ) 1 A<r- 1 + ... + ^ m) 

=aA = , /c = 1, g , 



(A - rj)(A - ri)yfR 2 g +1 (T,\) 



which are identical and imply af 1 \r) = bj(r), j — 0, 1..., g — 1. Therefore the ratio 



a, 



(m) 



r 



^ (r,A) -aa<™> J a,aM " Jfo (r) ; (22) 



does not depend on A and to check the fulfilment of the equality (|IzD for the functions 
Wj(X; r) defined by (|15D, (0) it is suffice to check it for any of terms of the expansion 
flTEl). The equation is obviously satisfied by the second term which gives the simplest 
nontrivial solution Wj = Vj. Hence the function (pl|) really is the generating function 
for the frequences of the KdV hierarchy. As a consequence, the generating equation 
for the phase wave number conservation laws of the KdV hierarchy has the form 

d T k^ + d x ^-A 3 / 2 =0. (23) 

3 General solution of the hodograph equations 

Now we are ready to construct the general solution to the hodograph equations (|T2|). 
Remind that it has to contain 2g + 1 arbitrary functions of one variable. So using the 
generating function Wj(r, A) as a kernel of convolution we arrive at the general solution 
in the form 

29+1 „ 

W j (r)= £ i^(r;%(A)dA, (24) 



k=1 k 

where 0fc(A) and are arbitrary functions and contours respectively on the Riemann 
surface T. The important restriction that has to be imposed upon 0fc's and A^s 
follows from the pure hyperbolicity of the modulational system (|8|) [9] and the form 
of the hodograph solution ([0]), namely, all W/s must be real. We remark that this 
requirement is not contained in the system (^) itself. Moreover, the contours A^ have 
to lie off the branch points. Finally, using the "potential" representation (T5[) , (|19"D 
we arrive at the solution 

W A*) = > j = 1,2,. ..,2^ + 1, m = l,...,g, (25) 



djk {m \r) 



where 



29+1 J 

Co (m) (r)= Uk(X)dnW . (26) 



Here "new" arbitrary 0^'s correspond to A 3,/2 0fc(A) from fl2~4]). Arbitrary constants 
0^ are real therefore the integration is accomplished along the axes in real r-space 
(this corresponds to a special choice of A^s which can be tightened to a real axis on 
the Riemann surface). Choosing </>fc(A) , one can obtain all the commuting flows 
described in [21,23]. 

Now we shall discuss briefly connection between the solution (p5|) and the Euler- 
Darboux-Poisson type systems which were investigated much earlier by Eisenhart [24] 
for the case of three independent variables (g = 1) and have arisen recently when 
studying the finite-gap averaged KdV equation [25,14,15,16,10]. In fact, the solution 
can be represented in the form 



UJ 



r) = E4 mj W/,(r) (27) 
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where ae»(r) are the normalization coefficients defined by (|7|) and 

2 fi ir / A«-y fc (A) . 

M r ) = E / / = rfA » * = (2? 



are the Eisenhart type solutions of the overdetermined compatible system of the Euler- 
Darboux-Poisson equations 

2(r i -r i )dJ/ <);f-<),f. i^j, i,j = l,...,2g + l. (29) 

The solution (fffi) generally is unlimited when any two of neighboring invariants r 3 - 
coalescing. For the further consideration we shall need the solution which is limited in 
all the domain of definition r\ < ... < r^g+i- The requirement of boundedness reduces 
the number of arbitrary functions in the general solution. We first demonstrate the way 
to construct the limited solution for the simplest single-phase case (see Refs. [14,15]). 
We start from the representation (pip which takes in the single-phase case the form 

3 Tk 

(30) 



k ~ l a k 



r-2 



dn = , k d\ , / (in = i . 



>R 3 (r, A) 

To provide boundedness of (50) at the singular surfaces r 2 = rj and r 2 = r 3 one 



splits the domain of definition D: r\ < r 2 < r 3 into two parts D\ : £ < r 2 < r 3 and 
D2 : ri < r 2 < £ with the deleted point (surface) r 2 = £, where ri < £ < r 3 is some 
real parameter on the Riemann surface. Each of these subdomains contains only the 
one singularity. To obtain the limited solution one should put a\ = a 2 = a 3 = £ and 

2 = -0 3 for L>i , 

(31) 

02 = -0i for D 2 ■ 

By this means the resulting solution is characterized by two (instead of initial three) 
arbitrary functions 0i(A) and 03(A). Introducing new (real) arbitrary functions ipi(X) = 
—0i (A) , ip2 (A) = \/— 103(A) we present the general bounded continuous solution in the 
form 

I V>i(A)dft + v 73 ! J ^ 2 (A)dft , for r x < r 2 < £ 

w(r) = (32) 

£ ^3 

| ^(A)dfi + v 73 ! | ^(A)rffi , for £ < r 2 < r 3 



'1 
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or more compactly 



min(r 2 ,0 r 3 

w(r)= / ^ 1 (A)c/fi + v /z T / i) 2 (\)dVL . (33) 

The obtained solution can be singular only at two points which are outside the domains 
D\ and D 2 , namely, r 2 = r\ = £ and r 2 = r% = £. These singularities are removed by 
the additional requirement 

^i(O=^2(O=0. (34) 

We remark that the solution ([33]), ([34]) generally is nonanalytic at the plane r 2 = 
£. This, however, does not contradict to a hyperbolic nature of the problem under 
consideration. Such a weak singularity was first revealed by Gurevich and Pitaevskii 
[13] in the numerical solution for the problem of cubic-like breaking (the exact solution 
of this problem, however, is analytic and does not contain any singularities [26,15,27]). 
The solution ([33]), ( |3"4"D can be easily generalized to the multiphase case. The result 
is: 



r ) = W , m = l,-,g, 



k=l 



where (35) 

min(r 2k £k) »*2fc+i 

4 m V) = / ^ 2k -i(\)dn^ + V^T J ^ 2k (\)dn^ , 

^2fe-l max(r 2 k,£k) 



^2fe-l(6) = i>2k(£k) = , 

r 2 k-i < Cfe < r 2k+i ■ (36) 

One can see that this solution is continuous and limited in all the domain of definition 
and generally contains g weak singular points r 2 k — , k — 1, ■■-,g. We remark in 
addition that formula (|35l) represents in some sense a decomposition of the g-phase 
solution into g single-phase ones in the hodograph space. 

The analytic limited solutions which were obtained in [10] in the form of infinite 
series can be picked out from the general representation (|35| ) by imposing the additional 
restriction 

^ 2 fc-i(A) = ^2fc(A) = if)(\) , k = l,...,g, (37) 
where is a real analytic function with g real zeros 



4 Boundary value problem solution 

Now we formulate the boundary problem to the system ([12]) corresponding to the 
initial value problem (HD, @. To avoid geometrical complifications we shall consider 



7 



only monotonically decreasing initial data (|I]). Evolution of such a curve with the only 
breaking point ( u"(xb r ) = ) is described completely (in the weak limit sense [1,2,6]) 
by the single-phase averaged system (|8|) with the boundary conditions formulated 
first by Gurevich and Pitaevskii [13]. The more recent investigations [1,2,16,28] have 
shown that these boundary conditions arise naturally in the zero dispersion limit of 
the KdV equation and are in the continuity of the Riemann invariants rj on the phase 
transition boundaries where r 2 k = r 2 k-i {j ^ 2k, 2k — 1) or r 2 k = r 2 k+i (j 7^ 2fc, 2k + 1). 
Earlier this problem was investigated by Avilov and Novikov [29] numerically in the 
framework of the theory of system of the hydrodynamic type. The effective analytic 
solutions of the Gurevich - Pitaevskii problem for all single-phase evolving (including 
nonanalytical and hump-like ) data were obtained in Refs. [15,30,31]. It is clear that 
every additional breaking introduces a new phase (gap) into consideration and therefore 
the genus of the problem g coincides with the number of breaking points on the initial 
curve [6]. By this means the number of dynamical variables in the problem is 2g + 1. 
At some moment t over i ap the evolving oscillatory phases will begin to overlap forming 
the multiphase solution region. To be specific we consider for the beginning the case 
when the initial curve Uq{x) has two distinct breaking points X\ and x 2 (Fig. 2a): 

u' Q (x) < , u'q(xi) = , u'q(x 2 ) = , 

v%(x x ) > , <'(z 2 ) > , (38) 
u (xi) = Ci , u (x 2 ) = c 2 . 

We assume in addition that the initial data are such that the breakings appear distinctly 
in space (i.e. that there are not points like A in Fig. 1). This assumption is made only 
for convenience of consideration and will be omitted afterwards. To obtain boundary 
conditions to the system (|12"D corresponding to the initial data (P, (|38D we apply 
to the one-phase stage of the evolution when two single-phase modes do not overlap 
yet (Fig. 2b). At this stage the global solution of (|12"D , which is generally defined 
in five (2g + 1) - dimensional r-space, lies on two (g) separated three-dimensional 
sections (r l5 r 2 , r 3 ), {r 2g _i, r 2g , r 2g+ i) (Fig. 3) (see Ref.[15]) with common diagonal 
surfaces r 2 = r 3 = r 4 , ...,r 2g _ 2 = r 2g _i = r 2g . These diagonal surfaces are the domains 
of definition of the zero-phase solution which separates the one-phase domains at the 
stage considered. Due to hyperbolicity of the hodograph equations (which is inherited 
from the hyperbolicity of the KdV-Whitham system) these 3-D domains do not interact 
and can be considered independently. We remark that the initial data and their smooth 
evolution until the first breaking are defined, as it will be shown below, on the pure 
diagonal r 1 = r 2 = ... = r 2g+ i in hodograph space [16,10]. 

Now we recall the important general fact proven in [10], which provides existence 
of the global solution to the initial value problem for the KdV-Whitham system: the 
uniform commuting flows Wj(r) (including Vj(r), of course) do not depend on multiple 
invariants r 2 k = r 2 k+\, or r 2 k = r 2 k-\ if j ^ 2k, 2k + 1 or 2k — 1, and, moreover, they 



8 



turn into their analogs for the g — 1 genus case: 



Wi 



{l"l,-;T 2 fc_l, P, (3, T 2k+2, 2g+l) 



(39) 



W 4 
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- 1 



r2jfe-i,r 2 jfc + 2, ...,r 2g +i), j^k,k + l 



and the analogous assertion is true for r 2k = r 2 fc_i, j ^ 2k, 2k — 1. As the obvious con- 
sequence of the equations (O) these relationships take place for all (not only uniform) 



Wj^s. This fact allows us to formulate boundary conditions at the phase transition 
surfaces r 2k = r 2k -i and r 2k = r 2 k+i, k = 1, ...,g which are single-phase - zero-phase 
transitions at the considered stage of evolution. Comparison of the hodograph solution 
( ]TT| ) with the solution of the Hopf equation (zero-phase averaged KdV) 

x - Qrt = W{r) , (40) 

where W(r) is the inverse function to the initial curve uo(x) ([[]), Q38|) (see Fig. 2a) 
shows that odd branches of the multivalued curve rj(x,r), j = 1, ...,5; r\ < ... < r$, 
have to be matched with the appropriate branches of the inverse function W(r) (see 
Fig. 2b). We shall demonstrate in detailes how to accomplish this matching at the 
boundaries where r 2 = r 3 and r 2 = r\ (the right-hand single-phase wave in Fig. 2b). 
The rest matching can be accomplished in an analogous way. 

As a first step we make a passage to the single-phase regime in the two-phase 
desired solution (|I2"D with certain Wj(ri, r 5 ). This corresponds to passage to times 

tbr ^ t <i t overlap- 
Using (|39"D we have the transition 



Wi (r 1 ,r 2 ,r 3 ,C,C) =W j {r x ,r 2 ,r z ) , j = 1,2,3, (41) 
9 = 2 g = l 

where ( is an arbitrary real parameter. 

Now the boundary conditions for the function Wj(ri,r 2 ,r 3 ) can be obtained by 
means of the usual Gurevich-Pitaevskii matching (see Refs. [13-15]) for the single- 
phase case, namely, at the edges where r 2 = r 3 and r 2 = n we have, taking into 
account (^), the single-phase - zero-phase transitions which drive us to the solution 
of the Hopf equation fliOD with the initial data ([!]), fl38[) that is: 



when r 2 = r 3 W\ 



when r 2 = r\ W% 



(ri,C!,ci) 
9 = 1 



(ci,ci,r 3 ) 
9 = 1 



(n) 

9 = 



(rs) 
^ = 



Win) 



(42) 



W(r 3 ) 
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Here following the ideology of [14,15] we have placed the coalesced invariants into 
the breaking point c\ which in a natural way marks the domain of definition of the 
sought solution with respect to the variables r\ and r 3 , namely (see Fig. 2b) 

n < ci < r 3 (43) 



Using arbitrarity of the parameter ( in (|4l| ) we present the boundary conditions for 
the desired W/s in the final form 

Wi(ri,ci,ci,ci,ci) = Win) , 

(44) 

W^(ci,ci,r 3 ,r3,r 3 ) = W(r 3 ) , 

An analogous consideration for the left-hand single-phase wave gives two more condi- 
tions 

W 3 (r 3 ,r 3 ,r 3 ,c 2 ,c 2 ) =W(r 3 ) , 

(45) 

W 5 (c 2 , c 2 , c 2 , ca, r 5 ) = W(r B ) , 

and, obviously, 

r 3 < c 2 < r 5 . (46) 



So we have four (2g) boundary conditions of the Goursat type for the system fll2|) 
which is, as it was mentioned above, specified by five (2g + 1) functions of one variable. 
The rest condition is the boundedness of the W/s in all the domain of definition which 
provides the boundedness of the hodograph solution (|12"D at finite x, t. The limited 
solution is given by formulas (p5]), ( j35|) and is parametrized precisely by 2g arbitrary 
functions of one variable. Substituting this solution into the boundary conditions (|44"1), 
fl45] ) one arrives after some calculations at four linear Abel equations resolving of which 
gives a simple representation for the unknown functions ^(A) in (35). 

^2fc-i(A) = — — / —==dx , 

A 



(47) 



«A) = ^/-^=^, fc = l,2. 
27r 7 v A — a; 



It should be noted that the arbitrary parameters in fl3"5|) have been chosen as follows 
(compare (||) with (||), (gBJ) ): 

6 = ci , 6 = c 2 (48) 



to satisfy the conditions (36). 
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Examination of the general case of evolution of the initial perturbation with g 
breaking points is completely analogous and does not involve any additional difficulties. 
The result is: the representation ( W7h is valid for any g. Thus, the arbitrary functions 
■0fc(A) in the general limited solution (^) are the Abel transformations of the parts of 
the inverse to the initial perturbation function W(r) which are divided by the breaking 
points. We emphasize that this result is valid only for monotonic initial data which have 
the one-valued inverse function. Evolution of the hump-like or large-scale oscillating 
initial perturbation requires an additional consideration which will be accomplished 
in the following papers. Appropriate results for the single-phase case can be found in 
Refs. [30,31]. 

So we have found the solution describing the motion of the multivalued curve rj (x, t) 
proceeding from the single-phase matching. Let us examine the initial value problem 
for the Whitham system (8) corresponding to the solved boundary problem. Making 
the consequent phase transitions like ([|l|) in the hodograph solution ([ll]) we arrive 
putting the arbitrary parameter ( equal to rj at the system 

•'• \'A'\j '•,)/ U'(0 0) (49) 

which is the system of the implicit solutions of 2g + 1 Hopf equations (zero-phase 
averaged KdV ) with the same initial data 

x-6r j t = W{r j ) j = 1, 2g + 1. (50) 

Therefore, both initial data and their smooth evolution are given on the pure diagonal 
r i = r 2 — ■■■ — r 2g+i in the real r-space. By this means the desired initial value problem 
for the Whitham-KdV system has the form 

ri(x, 0) = r 2 (x, 0) = ... = r 2g+1 (x, 0) = u (x) . (51) 

It should be noted that our consideration really did not require any assumptions of the 
distinct breakings so the latter can be omitted. 

Now we determine the motion of the phase transition boundaries which are the 
multiple caustics for the obtained solution. Really, it follows from the potential repre- 
sentation (H) that V 2 k(r) coalesces with either V2fc_i(r) or V 2 fe + i(r) when r 2 k coalesces 
with either r 2 fc_i or r 2 fc+i respectively. By this means the lines Xj(t) of the phase tran- 
sition in the x, t-plane have the multiple characteristic direction in every point and are 
given by the ordinary equations 

37 = Vj multipleijc) (52) 

which should be considered together with the hodograph solution (JTTJ) , (£25|), (ft5]), 

Let us discuss briefly the evolution of the Riemann surface topology on the obtained 
family of solutions. The Riemann surface of genus g is topologically equivalent to a 
sphere with g handles. These handles evolve both in space and in time and some of 
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them degenerate into poles on the phase transition boundaries in (x, t)-plane. Tending 
time to zero we arrive at the initial Riemann surface which represents a sphere with 
g singularities (poles). By this means we conclude that the topology of the evolving 
Riemann surface from the very beginning is given not only by the KdV equation itself 
but by the whole Cauchy problem formulation. 



5 The fundamental solution to the initial value prob- 
lem. 

We give another representation for the solution to the initial value problem for the 
Whitham-KdV system. Substituting the Abel transformations (H7J) into the solution 
( 35 ) and interchanging the order of integration we find 



T2g + 1 

u im) = j W(x)K^ m \v,x)dx (53) 
n 

where K^ m \r,x) is the fundamental solution (Riemann kernel function [32] ) to the 
initial value problem (g), (pT|) : 



min(rok.x) , \ 

(-1) / - for r 2fe _i < x < c k , 

J vx — A 

27rir (m) (r,x) = (54) 

' -l) k j - for c k < x < r 2k +i ■ 



Vx - X 



6 Summary and Conclusions 

The generating function of the uniform commuting flows for the multiphase averaged 
KdV equation (Whitham-KdV hierarchy) has been obtained in a direct way as de- 
pending on an additional parameter A solution of the linear generalized hodograph 
equations. It has a potential form and the potential (generalized frequency) turns out 
to be proportional to the coefficient of the holomorphic basis differential on the Rie- 
mann surface of genus g, where g coincides with the number of phases. The general 
solution in hodograph space is constructed as a superposition of 2g + 1 contour inte- 
grals of this generating function multiplied by arbitrary functions of one variable. The 
physical solutions corresponding to the Cauchy problem with monotonically decreas- 
ing initial data with g breaking points for the initial KdV equation are distinguished 
by the requirement of boundedness and depend on 2g arbitrary functions. Resolving 
of the appropriate boundary problem to the Whitham equations in hodograph space 
shows that these functions are the linear Abel transformations of the parts of the initial 
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curve divided by breaking points. Evolution of the Riemann surface topology on the 
obtained family of solutions is discussed and it is shown that "initial" Riemann sur- 
face represents a sphere with g singularities (poles). The initial value problem for the 
Whitham-KdV equations has been formulated and its fundamental solution (Riemann 
kernel function) has been constructed. 
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Figure Captions 

• 1. Splitting of the (x, £)-plane in the zero dispersion limit. 

• 2. Riemann invariants behaviour for the case g — 2 at a) t — b) £ br < £ < 

^overlap 

• 3. Domain of definition of the single-phase mode: r 2 k-i < ?"2fc < ^2fc+i 
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